Nifuroxazide suppresses PD-L1 expression and enhances the efficacy of radiotherapy in hepatocellular carcinoma.

Radiation therapy is a primary treatment for hepatocellular carcinoma (HCC), but its effectiveness can be diminished by various factors. The over-expression of PD-L1 has been identified as a critical reason for radiotherapy resistance. Previous studies have demonstrated that nifuroxazide exerts antitumor activity by damaging the Stat3 pathway, but its efficacy against PD-L1 has remained unclear. In this study, we investigated whether nifuroxazide could enhance the efficacy of radiotherapy in HCC by reducing PD-L1 expression. Our results showed that nifuroxazide significantly increased the sensitivity of tumor cells to radiation therapy by inhibiting cell proliferation and migration while increasing apoptosis in vitro. Additionally, nifuroxazide attenuated the up-regulation of PD-L1 expression induced by irradiation, which may be associated with increased degradation of PD-L1 through the ubiquitination-proteasome pathway. Furthermore, nifuroxazide greatly enhanced the efficacy of radiation therapy in H22-bearing mice by inhibiting tumor growth, improving survival, boosting the activation of T lymphocytes, and decelerating the ratios of Treg cells in spleens. Importantly, nifuroxazide limited the increased expression of PD-L1 in tumor tissues induced by radiation therapy. This study confirms, for the first time, that nifuroxazide can augment PD-L1 degradation to improve the efficacy of radiation therapy in HCC-bearing mice.

FluPMT: Prediction of Predominant Strains of Influenza A Viruses Via Multi-task Learning.

Seasonal influenza vaccines play a crucial role in saving numerous lives annually. However, the constant evolution of the influenza A virus necessitates frequent vaccine updates to ensure its ongoing effectiveness. The decision to develop a new vaccine strain is generally based on the assessment of the current predominant strains. Nevertheless, the process of vaccine production and distribution is very time-consuming, leaving a window for the emergence of new variants that could decrease vaccine effectiveness, so predictions of influenza A virus evolution can inform vaccine evaluation and selection. Hence, we present FluPMT, a novel sequence prediction model that applies an encoder-decoder architecture to predict the hemagglutinin (HA) protein sequence of the upcoming season's predominant strain by capturing the patterns of evolution of influenza A viruses. Specifically, we employ time series to model the evolution of influenza A viruses, and utilize attention mechanisms to explore dependencies among residues of sequences. Additionally, antigenic distance prediction based on graph network representation learning is incorporated into the sequence prediction as an auxiliary task through a multi-task learning framework. Experimental results on two influenza datasets highlight the exceptional predictive performance of FluPMT, offering valuable insights into virus evolutionary dynamics, as well as vaccine evaluation and production.

Attenuated Salmonella carrying siRNA-PD-L1 and radiation combinatorial therapy induces tumor regression on HCC through T cell-mediated immuno-enhancement.

Hepatocellular carcinoma (HCC), the most prevalent type of aggressive liver cancer, accounts for the majority of liver cancer diagnoses and fatalities. Despite recent advancements in HCC treatment, it remains one of the deadliest cancers. Radiation therapy (RT) is among the locoregional therapy modalities employed to treat unresectable or medically inoperable HCC. However, radioresistance poses a significant challenge. It has been demonstrated that RT induced the upregulation of programmed death ligand 1 (PD-L1) on tumor cells, which may affect response to PD-1-based immunotherapy, providing a rationale for combining PD-1/PD-L1 inhibitors with radiation. Here, we utilized attenuated Salmonella as a carrier to explore whether attenuated Salmonella carrying siRNA-PD-L1 could effectively enhance the antitumor effect of radiotherapy on HCC-bearing mice. Our results showed that a combination of siRNA-PD-L1 and radiotherapy had a synergistic antitumor effect by inhibiting the expression of PD-L1 induced by radiation therapy. Mechanistic insights indicated that the combination treatment significantly suppressed tumor cell proliferation, promoted cell apoptosis, and stimulated immune cell infiltration and activation in tumor tissues. Additionally, the combination treatment increased the ratios of CD4+ T, CD8+ T, and NK cells from the spleen in tumor-bearing mice. This study presents a novel therapeutic strategy for HCC treatment, especially for patients with RT resistance.

Prediction of Antigenic Distance in Influenza A Using Attribute Network Embedding.

Owing to the rapid changes in the antigenicity of influenza viruses, it is difficult for humans to obtain lasting immunity through antiviral therapy. Hence, tracking the dynamic changes in the antigenicity of influenza viruses can provide a basis for vaccines and drug treatments to cope with the spread of influenza viruses. In this paper, we developed a novel quantitative prediction method to predict the antigenic distance between virus strains using attribute network embedding techniques. An antigenic network is built to model and combine the genetic and antigenic characteristics of the influenza A virus H3N2, using the continuous distributed representation of the virus strain protein sequence (ProtVec) as a node attribute and the antigenic distance between virus strains as an edge weight. The results show a strong positive correlation between supplementing genetic features and antigenic distance prediction accuracy. Further analysis indicates that our prediction model can comprehensively and accurately track the differences in antigenic distances between vaccines and influenza virus strains, and it outperforms existing methods in predicting antigenic distances between strains.

Identification of and Mechanistic Insights into SARS-CoV-2 Main Protease Non-Covalent Inhibitors: An In-Silico Study.

The indispensable role of the SARS-CoV-2 main protease (Mpro) in the viral replication cycle and its dissimilarity to human proteases make Mpro a promising drug target. In order to identify the non-covalent Mpro inhibitors, we performed a comprehensive study using a combined computational strategy. We first screened the ZINC purchasable compound database using the pharmacophore model generated from the reference crystal structure of Mpro complexed with the inhibitor ML188. The hit compounds were then filtered by molecular docking and predicted parameters of drug-likeness and pharmacokinetics. The final molecular dynamics (MD) simulations identified three effective candidate inhibitors (ECIs) capable of maintaining binding within the substrate-binding cavity of Mpro. We further performed comparative analyses of the reference and effective complexes in terms of dynamics, thermodynamics, binding free energy (BFE), and interaction energies and modes. The results reveal that, when compared to the inter-molecular electrostatic forces/interactions, the inter-molecular van der Waals (vdW) forces/interactions are far more important in maintaining the association and determining the high affinity. Given the un-favorable effects of the inter-molecular electrostatic interactions-association destabilization by the competitive hydrogen bond (HB) interactions and the reduced binding affinity arising from the un-compensable increase in the electrostatic desolvation penalty-we suggest that enhancing the inter-molecular vdW interactions while avoiding introducing the deeply buried HBs may be a promising strategy in future inhibitor optimization.

Mechanistic Origin of Different Binding Affinities of SARS-CoV and SARS-CoV-2 Spike RBDs to Human ACE2.

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein (RBDCoV2) has a higher binding affinity to the human receptor angiotensin-converting enzyme 2 (ACE2) than the SARS-CoV RBD (RBDCoV). Here, we performed molecular dynamics (MD) simulations, binding free energy (BFE) calculations, and interface residue contact network (IRCN) analysis to explore the mechanistic origin of different ACE2-binding affinities of the two RBDs. The results demonstrate that, when compared to the RBDCoV2-ACE2 complex, RBDCoV-ACE2 features enhanced dynamicsand inter-protein positional movements and increased conformational entropy and conformational diversity. Although the inter-protein electrostatic attractive interactions are the primary determinant for the high ACE2-binding affinities of both RBDs, the significantly enhanced electrostatic attractive interactions between ACE2 and RBDCoV2 determine the higher ACE2-binding affinity of RBDCoV2 than of RBDCoV. Comprehensive comparative analyses of the residue BFE components and IRCNs between the two complexes reveal that it is the residue changes at the RBD interface that lead to the overall stronger inter-protein electrostatic attractive force in RBDCoV2-ACE2, which not only tightens the interface packing and suppresses the dynamics of RBDCoV2-ACE2, but also enhances the ACE2-binding affinity of RBDCoV2. Since the RBD residue changes involving gain/loss of the positively/negatively charged residues can greatly enhance the binding affinity, special attention should be paid to the SARS-CoV-2 variants carrying such mutations, particularly those near or at the binding interfaces with the potential to form hydrogen bonds and/or salt bridges with ACE2.

A Deep Learning Approach for Predicting Antigenic Variation of Influenza A H3N2.

Modeling antigenic variation in influenza (flu) virus A H3N2 using amino acid sequences is a promising approach for improving the prediction accuracy of immune efficacy of vaccines and increasing the efficiency of vaccine screening. Antigenic drift and antigenic jump/shift, which arise from the accumulation of mutations with small or moderate effects and from a major, abrupt change with large effects on the surface antigen hemagglutinin (HA), respectively, are two types of antigenic variation that facilitate immune evasion of flu virus A and make it challenging to predict the antigenic properties of new viral strains. Despite considerable progress in modeling antigenic variation based on the amino acid sequences, few studies focus on the deep learning framework which could be most suitable to be applied to this task. Here, we propose a novel deep learning approach that incorporates a convolutional neural network (CNN) and bidirectional long-short-term memory (BLSTM) neural network to predict antigenic variation. In this approach, CNN extracts the complex local contexts of amino acids while the BLSTM neural network captures the long-distance sequence information. When compared to the existing methods, our deep learning approach achieves the overall highest prediction performance on the validation dataset, and more encouragingly, it achieves prediction agreements of 99.20% and 96.46% for the strains in the forthcoming year and in the next two years included in an existing set of chronological amino acid sequences, respectively. These results indicate that our deep learning approach is promising to be applied to antigenic variation prediction of flu virus A H3N2.

Exploring the Cold-Adaptation Mechanism of Serine Hydroxymethyltransferase by Comparative Molecular Dynamics Simulations.

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to µs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.

The Energetic Origin of Different Catalytic Activities in Temperature-Adapted Trypsins.

Psychrophilic enzymes were always observed to have higher catalytic activity (k cat) than their mesophilic homologs at room temperature, while the origin of this phenomenon remains obscure. Here, we used two different temperature-adapted trypsins, the psychrophilic Atlantic cod trypsin (ACT) and the mesophilic bovine trypsin (BT), as a model system to explore the energetic origin of their different catalytic activities using computational methods. The results reproduce the characteristic changing trends in the activation free energy, activation enthalpy, and activation entropy between the psychrophilic and mesophilic enzymes, where, in particular, the slightly decreased activation free energy of ACT is determined by its considerably reduced activation enthalpy rather than by its more negative activation entropy compared to BT. The calculated electrostatic contributions to the solvation free energies in the reactant state/ground sate (RS/GS) and transition state (TS) show that, going from BT to ACT, the TS stabilization has a predominant effect over the RS stabilization on lowering the activation enthalpy of ACT. Comparison between the solvation energy components reveals a more optimized electrostatic preorganization to the TS in ACT, which provides a larger stabilization to the TS through reducing the reorganization energy, thus resulting in the lower activation enthalpy and hence lower activation free energy of ACT. Thus, it can be concluded that it is the difference in the protein electrostatic environment, and hence its different stabilizing effects on the TS, that brings about the different catalytic activities of different temperature-adapted trypsins.

Molecular dynamics simulations reveal distinct differences in conformational dynamics and thermodynamics between the unliganded and CD4-bound states of HIV-1 gp120.

The entry of human immunodeficiency virus type I (HIV-1) into host cells is initiated by binding to the cell-surface receptor CD4, which induces a conformational transition of the envelope (Env) glycoprotein gp120 from the closed, unliganded state to the open, CD4-bound state. Despite many available structures in these two states, detailed aspects on the dynamics and thermodynamics of gp120 remain elusive. Here, we performed microsecond-scale (µs-scale) multiple-replica molecular dynamics (MD) simulations to explore the differences in the conformational dynamics, protein motions, and thermodynamics between the unliganded and CD4-bound/complexed forms of gp120. Comparative analyses of MD trajectories reveal that CD4 binding promotes the structural deviations/changes and conformational flexibility, loosens the structural packing, and complicates the molecular motions of gp120. Comparison of the constructed free energy landscapes (FELs) reveals that the CD4-complexed gp120 has more conformational substates, larger conformational entropy, and lower thermostability than the unliganded form. Therefore, the unliganded conformation represents a structurally and energetically stable "ground state" for the full-length gp120. The observed great increase in the mobility of V1/V2 and V3 along with their more versatile movement directions in the CD4-bound gp120 compared to the unliganded form suggests that their orientations with respect to each other and to the structural core determine the differences in the conformational dynamics and thermodynamics between the two gp120 forms. The results presented here provide a basis by which to better understand the functional and immunological properties of gp120 and, furthermore, to deploy appropriate strategies for the development of anti-HIV-1 drugs or vaccines.

Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases.

To investigate the role of electrostatics in different temperature adaptations, we performed a comparative study on subtilisin-like serine proteases from psychrophilic Vibrio sp. PA-44 (VPR), mesophilic Engyodontium album (Tritirachium album) (PRK), and thermophilic Thermus aquaticus (AQN) using multiple-replica molecular dynamics (MD) simulations combined with continuum electrostatics calculations. The results reveal that although salt bridges are not a crucial factor in determining the overall thermostability of these three proteases, they on average provide the greatest, moderate, and least electrostatic stabilization to AQN, PRK, and VPR, respectively, at the respective organism growth temperatures. Most salt bridges in AQN are effectively stabilizing and thus contribute to maintaining the overall structural stability at 343 K, while nearly half of the salt bridges in VPR interconvert between being stabilizing and being destabilizing, likely aiding in enhancing the local conformational flexibility at 283 K. The individual salt bridges, salt-bridge networks, and calcium ions contribute differentially to local stability and flexibility of these three enzyme structures, depending on their spatial distributions and electrostatic strengths. The shared negatively charged surface potential at the active center of the three enzymes may provide the active-center flexibility necessary for nucleophilic attack and proton transfer. The differences in distributions of the electro-negative, electro-positive, and electro-neutral potentials, particularly over the back surfaces of the three proteases, may modulate/affect not only protein solubility and thermostability but also structural stability and flexibility/rigidity. These results demonstrate that electrostatics contributes to both heat and cold adaptation of subtilisin-like serine proteases through fine-tuning, either globally or locally, the structural stability and conformational flexibility/rigidity, thus providing a foundation for further engineering and mutagenesis studies.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations.

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.

Insights into Protein-Ligand Interactions: Mechanisms, Models, and Methods.

Molecular recognition, which is the process of biological macromolecules interacting with each other or various small molecules with a high specificity and affinity to form a specific complex, constitutes the basis of all processes in living organisms. Proteins, an important class of biological macromolecules, realize their functions through binding to themselves or other molecules. A detailed understanding of the protein-ligand interactions is therefore central to understanding biology at the molecular level. Moreover, knowledge of the mechanisms responsible for the protein-ligand recognition and binding will also facilitate the discovery, design, and development of drugs. In the present review, first, the physicochemical mechanisms underlying protein-ligand binding, including the binding kinetics, thermodynamic concepts and relationships, and binding driving forces, are introduced and rationalized. Next, three currently existing protein-ligand binding models--the "lock-and-key", "induced fit", and "conformational selection"--are described and their underlying thermodynamic mechanisms are discussed. Finally, the methods available for investigating protein-ligand binding affinity, including experimental and theoretical/computational approaches, are introduced, and their advantages, disadvantages, and challenges are discussed.

Comparison of bimanual and micro-coaxial phacoemulsification with torsional ultrasound.

PURPOSE: To evaluate the efficacy and postoperative outcomes of bimanual and micro-coaxial phacoemulsification with torsional ultrasound. METHODS: In this prospective randomized series, eyes with age-related cataract were randomly divided into bimanual phacoemulsification and micro-coaxial phacoemulsification groups. Data included preoperative and postoperative corrected distance visual acuity (CDVA), central/temporal corneal thickness, endothelial cell density and intraoperative ultrasound time, cumulative dissipated energy and balanced salt solution volume. RESULTS: The study evaluated 89 patients (89 eyes). When compared between micro-coaxial and bimanual phaco groups, ultrasound time (57 ± 24 second versus 85 ± 40 second, p < 0.01), cumulative dissipated energy (10 ± 7 versus 14 ± 7, p = 0.01) and balanced salt solution volume (55 ± 19 ml l versus 75 ± 20 ml, p < 0.01) were significantly lower in micro-coaxial phaco group. At 1 day, there were no statistically significant differences in CDVA (p = 0.68) or central corneal thickness (p = 0.48) between two groups; temporal corneal thickness was statistically significantly thicker in bimanual phaco group (1247 ± 123 µm) than that in micro-coaxial phaco group (1108 ± 131 µm, p = 0.01). There were no statistically significant differences in CDVA, central/temporal corneal thickness, or endothelial cell density between two groups 30 days postoperatively (all p values >0.05). CONCLUSION: With rapid visual rehabilitation and comparable endothelial cell loss, both bimanual and micro-coaxial phacoemulsification procedures with torsional ultrasound are efficient for moderate nuclear cataract extraction.

Torsional ultrasound efficiency under different vacuum levels in different degrees of nuclear cataract.

PURPOSE: To evaluate the efficiency of torsional ultrasound (US) under different vacuum levels in eyes with cataract with different degrees of nuclear density. SETTING: Cataract Service, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China. METHODS: In this prospective randomized series, eyes with age-related cataract were divided into 3 groups based on nuclear density. Eyes in each group were randomly assigned to 1 of 3 subgroups and treated with phacoemulsification under vacuum of 250 mm Hg (subgroup 1), 450 mm Hg (subgroup 2), or 600 mm Hg (subgroup 3). Data included preoperative and postoperative corrected distance visual acuity (CDVA), endothelial cell density, and central corneal thickness (CCT) and intraoperative US time and cumulative dissipated energy (CDE). RESULTS: The study evaluated 998 patients (1073 eyes). In each group, there was a statistically significant decrease in US time, CDE, and postoperative endothelial cell loss with increased vacuum (P<.05). At 1 day and 7 days, the CDVA was statistically significantly better in subgroups 2 and 3 (P<.05) and the CCT was thinner when vacuum was higher. There was no statistically significant difference in CDVA and CCT between subgroups 30 days postoperatively (P>.05). CONCLUSIONS: Torsional US with a high vacuum level was safe for cataract extraction. With less US energy and endothelial cell loss, torsional US was more efficient than with higher vacuum levels with lower levels.

Early changes in clear cornea incision after phacoemulsification: an anterior segment optical coherence tomography study.

PURPOSE: We aimed to investigate the ultrastructures of clear corneal incisions 24 hours after phacoemulsification using anterior segment optical coherence tomography (AS-OCT). METHODS: Sixty eyes of 60 patients scheduled for cataract surgery were randomly selected. All eyes underwent the same phacoemulsification procedure carried out by one experienced surgeon. Two-plane temporal clear corneal tunnel incisions were performed. The same types of intraocular lens and implant system were used in all patients. Images of the ultrastructures of the corneal incisions were taken using AS-OCT before and 1 day after surgery. RESULTS: Corneal thickness increased at the incision site after surgery in all eyes (p < 0.001). Epithelial bulla in the incision region was seen in two eyes (3%). Gaping at the internal aspect of the corneal wound was seen in 42 eyes (70%). Eyes with gaping had thicker localized cornea (p = 0.002). Descemet's membrane detachment was seen in 49 eyes (82%); this seemed to be associated with lower preoperative intraocular pressure (p = 0.01). CONCLUSIONS: Anterior segment OCT provides sensitive and detailed measurements of the ultrastructures in clear corneal incision.

Comparison of lens thickness measurements using the anterior segment optical coherence tomography and A-scan ultrasonography.

PURPOSE: To assess lens thickness (LT) measurements with anterior segment optical coherence tomography (AS-OCT) in comparison with A-scan ultrasonography (A-scan US). METHODS: Enrolled were 66 eyes of 48 phakic elderly volunteers aged > or = 50 years and 56 eyes of 56 young participants aged 18 to 40 years. LT was measured with the internal manual caliper tools in AS-OCT. The A-scan US measurements were based on an average of 10 consecutive automatic measurements. Reproducibility was assessed by three measurements each with AS-OCT and A-scan US independently obtained by two observers who were masked to one another's results. RESULTS: The failure rates of AS-OCT and A-scan US were 9.1% and 7.6%, respectively, in elderly subjects, but no failure was observed in young subjects. The LT values measured by AS-OCT were significantly greater than A-scan US; the paired difference was 0.135 mm in elderly and 0.101 mm in young subjects (P < 0.001). These differences did not correlate with the nuclear cataract grades (r = 0.078, P = 0.558). Intraobserver agreement on AS-OCT (95% limits of agreement [LoA]: -0.049 to +0.045 mm; ICC: 0.999) was better than A-scan US (95%LoA: -0.194 to +0.218 mm; ICC: 0.974). The 95% LoA of interobserver agreement using AS-OCT and A-scan were -0.084 to +0.073 and -0.278 to +0.239 mm, respectively, and the ICCs were 0.996 and 0.960, respectively. CONCLUSIONS: AS-OCT can be used to measure lens thickness in most eyes with clear or opacified lenses. It appears to be an alternative means of measuring lens thickness, particularly when a noncontact method is needed.